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Abstract -We present measurements on hybrid single-electron turnstiles with superconducting 
leads contacting a normal island (SINIS). We observe Andreev tunneling of electrons influencing 
the current plateau characteristics of the turnstiles under radio-frequency pumping. The data is 
well accounted for by numerical simulations. We verify the dependence of the Andreev tunneling 
rate on the turnstile's charging energy. Increasing the charging energy effectively suppresses the 
Andreev current. 



Introduction. — At the present time, there does not 
exist a quantum current standard though such standards 
exist for both voltage and resistance. A strong candidate 
for such a standard is a turnstile with a small normal metal 
island connected by tunnel junctions to superconducting 
leads (SINIS) [l]. These hybrid turnstiles pump electrons 
one at a time producing a well-defined current of / = e/ 
where e is the electron charge and / is the frequency of 
pumping. Understanding and eliminating error processes 
in these turnstiles is vital for realizing a quantum metro- 
logical triangle (QMT), a key goal in metrology [2]. In 
closing the QMT, the three standards of voltage, resis- 
tance and current would be compared against each other 
via Ohm's law. Closing the QMT will allow the most ac- 
curate comparison of the Josephson constant Kj = 2e/h 
and the von Klitzing constant Rk = h/e^ (i.e. the charge 
of the electron and Planck's constant) to date. 

There exist several potential quantum current standard 
candidates. The NIST seven junction pump has demon- 
strated a current accuracy of 1.5 parts in 10® [3j|4] but 
is limited to maximum currents of approximately 1 pA. 
Other candidates with the potential for metrological ac- 
curacy at metrologically relevant currents include semi- 
conductor tunable barrier pumps, charge-coupled device 
pumps |5], tunnel junction pumps |6}|8|, quantum dot 



pumps [9] , surface acoustic wave pumps 10 and quantum 
phase slip nanowire pumps 11 . In general, the accuracy 



ble to the accuracy of the quantum Hall and Josephson 
effect standards. A semiconducting quantum dot parallel 
pump recently demonstrated an accuracy of 1.5 in 10^ at 



54 pA 12 . The minimum required current magnitude for 
closing the metrological triangle is about 100 pA with an 
accuracy of about 1 part in 10®, which has not yet been 
achieved. 

The SINIS turnstile operates by using energy barriers 
to control the flow of electrons. The normal metal island 
is capacitively coupled to a gate electrode. By applying 
an appropriate voltage to the gate, single electrons can 
be added and removed from the normal metal island. In 
essence, the SINIS turnstile operates as a single-electron 
transistor, with the superconducting gap providing extra 
protection against unwanted tunneling of the electrons. 

There are multiple transport processes that can occur in 
a SINIS turnstile. The dominant one employed in charge 
pumping is sequential single electron tunneling through 
the insulating barrier. The dominant two-electron process 
causing errors is Andreev tunneling. In Andreev tunnel- 
ing, an electron in the normal metal is reflected as a hole 
(or a hole as an electron) at the interface of the NIS junc- 
tion. This forms (or removes) a Cooper pair in the super- 
conducting electrode and can be alternatively viewed as 
two electrons tunneling simultaneously across the insulat- 
It has been shown that errors arising from 



ing barrier 13 



of these quantum currents standards is not yet compara- 



sequential le-tunneling in turnstiles such as environmental 
activation can be effectively suppressed by proper flltering 
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including an on-chip capacitively coupled ground plane 



14 



Higher order, multiple-electron error processes not 
eliminated by this filtering may then be observed. Many 
experimental observations of Andreev currents have been 
reported previously 
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Andreev processes in devices 
with Coulomb blockade have been studied previously, both 
theoretically .IS, and experimentally 13 19 . Only re- 



cently were individual Andreev reflection events detected 



in NIS junctions 20 21 



Here wc demonstrate that the Andreev tunneling pro- 
cess is detectable in single-electron turnstile current pump- 
ing plateaus and that the resulting error effect can be min- 
imized by increasing the charging energy of the turnstile. 
The data is well accounted for by numerical simulations 
with reasonable parameter values and the observed behav- 
ior is physically intuitive. 
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Fig. 1: Experimental set-up (a) Layout for sample of type one. 
The leads are coupled to the ground plane by large ground 
plane couplers (GPC) yet the turnstile itself is off the ground 
plane, (b) A SEM micrograph of the first sample type showing 
the turnstile and combined RF gate. The GPC and ground 
plane are not shown. This is sample SI with junction area 50 
nm by 125 nm. (c) The second sample set-up with the samples 
in the patterned part of the ground plane. The device design 
is modified to include a separate RF and DC gate, (d) A SEM 
micrograph of the second sample type showing the turnstile 
and the DC gate on top of the RF gate. This is sample S4 
with junction area 60 nm by 70 nm. 



Experiment. — To form the SINIS turnstile, a nor- 
mal copper island is connected to two superconducting 
aluminum electrodes using double angle evaporation with 
an intermediate oxidation step. Two types of devices were 
fabricated using electron beam lithography. The first type 
of sample used one step of lithography and aligned the de- 
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Fig. 2: 10 MHz sine wave current pumping with simulations, 
(a) DC current vs voltage envelope measurement for sample SI 
(green points). The red and blue lines are simulations of the 
measurement, (b) Current pumping plateaus for low charging 
energy sample SI with a 10 MHz sine wave, (c) Close up of (b) 
showing the Andreev tunneling, (d) Close up of current pump- 
ing plateaus of sample S3 (high charging energy, low resistance) 
with a 10 MHz sine wave showing the absence of Andreev tun- 
neling due to higher charging energy. Bias voltages are: 160 
/iV (blue), 240 /iV (red) and 320 /iV (green) with a magenta 
dotted line at / = e/ and simulations shown as black dashed 
lines in figures b-d. 



vice with a ground plane as shown in fig. [TJa). Large pads 
ensured capacitive coupling between the device and the 
ground plane, strongly suppressing environmental noise 
14 20]. These large pads are marked ground plane cou- 



plers (GPC) and are electrically isolated from the ground 
plane by atomic layer deposition (ALD) grown aluminum 
oxide. This device design allows simple single sample fab- 
rication but is difficult to parallelize. A single gate line 
is used as both a DC gate and an RF gate. A scanning 
electron microscope (SEM) micrograph of the first type of 
sample is shown in fig. [l|^b). 

The second type of sample involves multiple steps of 
lithography as shown in fig. [l|c). In the first step, a 
ground plane is patterned to allow device fabrication on 
top of it. The turnstile sits in the patterned gap in the 
ground plane and only the leads are directly above it. Pat- 
terning the ground plane in this manner prevents forma- 
tion of short circuits when wire bonding. After lithogra- 
phy, an insulating layer of silicon oxide is deposited by 
plasma enhanced chemical vapor deposition (PECVD). 
A second step of lithography is performed to form the 
turnstiles and individual DC gates for each turnstile. The 
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Fig. 3: Close ups of frequency and waveform dependence of 
current pumping plateaus (a) Pumping sample SI with a 10 
MHz sine wave showing the enhanced tunneling caused by the 
Andreev process, (b) Pumping of sample S2 (high charging en- 
ergy, high resistance) with a 10 MHz sine wave showing the sup- 
pression of the enhanced Andreev tunneling. At high pumping 
amplitude, there is increased back tunneling compared to sam- 
ple S3 due to the higher resistance of the junctions, (c) Pump- 
ing sample SI with a 10 MHz square wave showing the different 
character of enhanced Andreev tunneling with a square wave 
drive, (d) Close up of current pumping plateaus of sample SI 
with a 50 MHz sine wave showing the smaller relative error at 
higher pumping frequency. The vertical scale has been changed 
compared to previous close-ups for clarity. The vertical scale 
in figures a-c are identical. The horizontal scale is not identical 
since the exact amplitude of RF drive applied to an individual 
device is not important. Bias voltages are: 180 /xV (blue), 220 
/iV (red) and 260 (green) with a magenta dotted line at 
I = ef and simulations shown as black dashed lines in all four 
figures. 



lithography used a three layer mask with a hard germa- 
nium layer which allowed smaller features than the two 
layer PMMA patterning used for the first type. The ad- 
vantage of this fabrication method is that the separate 
DC gates enables parallelization. Parallelization is essen- 
tial for metrological purposes as it is difficult to get large 

A SEM micro- 



enough currents from a single device 22 
graph of the second type of sample is shown in fig. [ijd). 

Measurements are done at a base temperature of ap- 
proximately 70 mK. A typical current vs voltage envelope 
measurement reveals the superconducting gap and single 
electron transistor behavior of the turnstile as shown in 
fig. [2|^a). To measure the pumped current of the turnstile, 
a fixed bias voltage is applied to one of the turnstile's bias 



leads while current is measured on the other lead. A sine 
wave of amplitude Vac (resulting in a normalized ampli- 
tude Ag = CgVAc/c) with a frequency / is applied to the 
RF gate. As Ag increases past a threshold value, the likeli- 
hood of transporting one electron through the turnstile in 
a single cycle approaches unity, giving rise to a quantized 
current plateau. 

We observe current quantization as shown in fig. [2]jb). 
Zooming in on the current plateau reveals characteristic 
excess current above the expected I = ef behavior due to 
Andreev tunneling. For samples with small charging en- 
ergy (Ec < A), as shown in fig. ^c), the current plateau 
is enhanced above the expected value particularly for low 
Ag values in the plateau. These deviations are suppressed 
for samples with large charging energies {Ec > A), as 
shown in fig. [2jd). Amplifier gain is corrected on the 10"'^ 
level to match the simulated plateaus. This does not influ- 
ence the interpretation of excess current on the 10~^ level. 
Qualitatively, as Ag is increased for the situation Ec < A, 
we encounter an energy threshold that permits Andreev 
tunneling before the single electron tunneling threshold is 
encountered, allowing enhanced current flow for low charg- 
ing energy samples. This charging energy dependence is 
indicative of the Andreev effect 
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Fig. [3] shows close ups of the current pumping plateaus 
for a different conditions. In fig. [sf^a), pumping of sample 
SI with a 10 MHz sine wave at a different range of biases 
than fig. [2] is shown. In fig. [sjjb) the pumping curves 
for the high charging energy, high resistance sample S2 
shows a reduction below the expected I = ef due to single 
electron back tunneling. This increased back tunneling at 
high pumping amplitude is a first order effect, which is 
due to the larger resistance of the high charging energy 
sample S2 compared to S3 [24]. Figs, [sjjc) and (d) show 
sample SI pumped with a 10 MHz square wave and a 50 
MHz sine wave. 

The Andreev effect is visible in a wide range of biases 
as shown in figs, and (c) andjsj^a) with no visible de- 
pendence on bias voltage. Although the exact position of 
the plateau depends on bias, there is no separation of the 
plateaus due to bias dependence indicating this is not an 
environmental activation effect 14 . By fabricating high 



charging energy samples with low resistance, we can limit 
both the Andreev and back tunneling effects in the turn- 
stiles as shown in figs. [SJd). This requires small junctions 
with highly transparent tunnel barriers i.e. junctions with 



low RC product as have been fabricated previously 24|25 

It should be noted that for device S3 with optimized 
tunnel barriers for suppressing Andreev and single electron 
back tunneling there is still significant structure visible at 
the 10"'^ level. The exact origin of this structure is not 
yet clear. It is possibly due to quasi-particle relaxation 
in the aluminum leads or noise still penetrating the fil- 
tering system. Preliminary simulations indicate that this 
structure can be accounted for if we consider overheated 
superconducting leads but future research is needed in this 
area. The focus here was to demonstrate control over the 
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Table 1: Sample Parameters. 



c) 2 



Name 




Rt (kr2) 


A (^cV) 


AcH (nm^) 


SI 


0.63 


160 


220 


30 


S2 


2.2 


1400 


210 


N/A 


S3 


1.4 


430 


210 


N/A 


S4 


0.75 


110 


210 


30 



Andreev error process. SINIS pumps are theoretically pre- 
dicted to achieve an accuracy of 1 part in 10^ with a cur- 
rent of 30 pA for a single pump 26 . The experimentally 
measured accuracy level to date is on the 10~^ level i.e. 
part in 10^. 



1 



Theory. — For the simulations shown in figs. [2] and 
[3] we calculate the average current flowing through the 
turnstile by numerically solving a master equation. The 
current through the left barrier (which is equal to the cur- 
rent through the right barrier in the steady state) is given 
by: 



(1) 



where P{'n, t) is the probability for finding n electrons on 
the central island at time t. The tunneling rate from 
the left lead to the island is TLi{n) = TsNiE'^i'n)) + 
2T ar{E'^^ (nj) and the tunneling rate from the island to 
the left lead is r/L(n) = rjvs(£^r(n)) + 2T ar{EJ;- (n)) . 
The factor of two in front of the Andreev rate Tar, 
compared to the single electron rates Fjvs.SAr, is to ac- 
count for the Andreev process transporting two electrons. 
The energy required to add or remove a single electron 
to/from the island (denoted by + or — ) from the left 
or right lead (denoted by L or R) is given by — 
±2Ec in-ng± 1/2) ± bVlm where Ec = e^/2C is the 
charging energy with C being the total capacitance of 
the island, n is the number of excess electrons on the 
island and Ug = CgVg/e is the normalized offset charge 
(the effective charge induced by applying a voltage Vg to 



the gate with capacitance Cg) 27 . The energy required 
to add or remove two electrons to/from the island (de- 
noted by -f-|- or ) in the Andreev process is given by 

= ±4Ec {n-ng±l)± 2eVLM- 

The time dependence of the probability is given by the 
master equation 



27 28 



dP[n,t) 
dt 



+r„_i,„P(n -l,t)+ r„+i,„P(n + 1, t) 
+r„_2.«P(n -2,t) + V,+2.nP{n + 2, t) 



(2) 
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Fig. 4: (a) Probability evolution during charge pumping. The 
probability P(n, t) is shown for different values of n: n = 
(green), n = 1 (blue), n = 2 (dotted black) and n = —\ (dot- 
ted red), (b) Rates of tunneling weighted by probabilities of 
occupation. The single electron tunneling rate through the left 
junction from n = to n = 1 is shown in solid green. The 
Andreev tunneling rate from n = to n = 2 is shown in dotted 
green and is on top of the relaxation rate from n = 2 to n = 1 
via single electron tunneling through the right electrode (thick 
black). Likewise, the single electron tunneling rate through the 
right junction from n = 1 to n = is shown in solid blue, the 
Andreev rate n = lton = — lis shown in dotted blue and is 
on top of the relaxation rate from n = — lton = through 
the left junction (thick red), (c) Stability diagram of sample 
SI. The blue diamonds are the thresholds for single electron 
tunneling, the red diamonds are the thresholds for Andreev 
tunneling and the thick black line is the pumping cycle used in 
(a) and (b). (d) Stability diagram of sample S2. The diamonds 
are reversed in order because of the higher charging energy of 
sample S2. 



where 



VsN{Et(n - 1)) + VsN{.El{n - 1)) 
V^s{El {n + !))+ TNsiEj,, (n + I)) 
TAR{E++{n~2))+TAR{E++{n-2)) (3) 
TARiEl-{n + 2)) + TARiE^'in + 2)) 

Equation[2]tracks the fiow of the probability. The higher 
order Andreev tunneling effects are included by consid- 
ering state changes from n to n ± 2 26 . The rate at 



Tn— l.n 
^n+2,n 



which single electrons tunnel from superconductor to nor- 
mal metal Tsn, or the rate at which single electrons tun- 
nel from the normal metal to the superconductor Fjvs, are 



obtained using first order perturbation theory 29 30 . 
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are given m equa- 
using second order perturbative 



The Andreev tunneling rates Tahi 
tion 3 in reference 
calculations. These rates depend on the charging energy 
Ec, the superconducting gap A, and the tunneling resis- 
tance Rt- Values for these parameters are obtained from 
DC current vs voltage envelope measurements (see inset 
of fig. [I]) which depend only on T^s^sn- ^ar has an ad- 
ditional fitting parameter since its magnitude is controlled 
by the quantity [fi/ Rxe^) /N where N is the effective num- 
ber of conduction channels, N = A/Ach [26]. A is the 
cross-sectional area of the junction estimated by scanning 
electron microscope (SEM) imaging (see fig. [ijfor relevant 
images and area estimates) and Aqh is the effective area 
of a conduction channel used as a fitting parameter. The- 
oretically, AcH ~ 2 nm^ though fitted values are typically 
much larger than this value and are interpreted as result- 
ing from inhomogeneities in the thickness of the oxide in 
the junctions [20]. In table [l] the fitting parameters Eq 
(in units of A), Rt, A and Ach are listed for each of 
the samples simulated. Note that Andreev rate fitting pa- 
rameters can not be determined for high charging energy 
samples S2 and S3 since the Andreev effect is suppressed 
below the measurement noise. 

The simulated probability evolution during charge 
pumping is shown in fig. |4|^a) for various charge states. 
The corresponding rates of tunneling weighted by the oc- 
cupation probabilities are shown in fig. Qb). These plots 
are for simulations of sample SI with Vl — Vr — 200 /iV, 
Ag = 0.74 and / = 10 MHz. Only forward tunneling is 
relevant here. In figs. |4|^c) and (d), we show the calculated 
stability diamonds for samples SI and S2 respectively. The 
minimal pair breaking energy for le-tunneling is A so the 
threshold for tunneling is E^ = A which is shown as the 
blue line. For the Andreev tunneling, this pair breaking is 
avoided so the threshold for tunneling is E^^ = and is 
shown by the red dashed line. 

Results and Discussion. — In figs. |4|^a) and (b), 
the probabilities and tunneling rates during one pumping 
cycle are shown for sample SI. As the transition rate from 
n = to n = 1 grows, the probability for being in state 
n = quickly drops in the beginning of the cycle and is 
replaced with a probability for being in state n — 1. At 
the same time, there is a detectable Andreev tunneling 
rate for going from state n = to n = 2 after which the 
n = 2 state quickly relaxes to n = 1 by single-electron 
tunneling. The low level occupation of the state n = 2 
is also detectable. Likewise, the opposite process where 
an electron leaves the turnstile is observable in the second 
part of the pumping cycle. As before, when the rate from 
n = 1 to n = grows, the most likely state becomes n = 0. 
The transition for n = 1 to 7i = — 1 by Andreev tunneling 
is also seen whereafter the n — —1 state quickly relaxes to 
n = state by single-electron tunneling. 

In figs. [2];b), [2];c) andjga), the pumping plateaus for 
six different biases and the corresponding fits are shown 
for low charging energy (Ec < A) sample SI. There is 



excellent agreement between the fits and the data. By 
looking at the stability diagram shown in fig. |4|^c) , we can 
qualitatively understand the observed effect. The stabil- 
ity diagram shows tunneling thresholds versus normalized 
bias voltage, V, and momentary gate charge, Ug. Starting 
in the diamond on the left, as the gate amplitude, Ag, is 
gradually increased, the Andreev threshold (shown as a 
red dotted line) is encountered first. Thus the Andreev 
process is most noticeable at low Ag. As Ag is increased 
further, we encounter the single electron tunneling thresh- 
old and this process quickly dominates, obscuring the An- 
dreev effect. This results in the enhanced current pumping 
plateau in figs. [2jc) and[3|^a). For comparison, as can be 
seen in fig. Qd), the high charging energy {Ec > A) sam- 
ple, S2, encounters the single electron tunneling threshold 
before the Andreev threshold. The Andreev effect is not 
observed as the electron has already tunneled before it en- 
ters the Andreev regime. Thus the plateau in figs, [sjd) 
and[2jd) is flat with no evidence of the enhanced tunneling 
effect. 

The lack of dependence on bias seen in figs. [2]Jc) and 
[sj^a) can similarly be explained by looking at the stability 
diagram. Changing the bias corresponds to changing the 
horizontal level of the black line in fig. ^c). However, 
the Andreev tunneling thresholds run parallel to the single 
electron tunneling thresholds so this has a very small effect 
on the observed current. 

In fig. [sjjc), the results from pumping with a 10 MHz 
square wave for sample SI are shown. The square wave 
is modeled with an exponential rise to the applied voltage 
with risetime 2.5 ns. Compared to the sine wave pumping 
of the same frequency, the enhanced current for the square 
wave is more peaked and dies away more quickly. At an 
amplitude of Ag — 0.75, the square wave pumped plateau 
is almost flat with a current less than 1.005e/ (i.e. closer 
to the ideal value of J = e/) while sine wave pumping 
with the same frequency and bias voltage has a plateau 
with significant slope at = 0.75 with a current greater 
than 1.005e/. This behavior can also be qualitatively un- 
derstood from the stability diagram. The square wave 
goes more immediately to the final Ag value thus spending 
less time in the vulnerable regions of the stability diagram. 
When we first encounter the Andreev threshold but before 
the single electron threshold is encountered, the Andreev 
process dominates and we see the sharp peak in current 
shown in fig. [sjjc). As Ag is increased, we encounter the 
single-electron tunneling process threshold and the escape 
process is dominated by it. There is only a small possibil- 
ity for the Andreev tunneling enhancement to occur since 
the square wave sweeps through that vulnerable region 
quickly. Thus the Andreev effect in the pumping plateau 
dies away more quickly for the square wave than for the 
sine wave. 

In figs. [Sjjd), the pumping plateaus for a 50 MHz sine 
wave are shown with the Andreev effect visible for low 
charging energy sample, SI. With higher frequency, more 
electrons are pumped per unit time producing a higher 
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current. Thus the relative accuracy improves compared 
to the lower frequency 10 MHz sine wave shown in fig. 
[3]^ a) since the absolute value of the excess Andreev current 
remains roughly the same. The Andreev tunneling contri- 
bution to the current remains roughly unchanged because 
the same process as described for the lower frequency sine 
wave takes place. This results in the more flat slope of the 

50 MHz pumping seen in figure [sj^d) but the Andreev pro- 
cess is still detectable. It should be noted that the same 
sample parameters listed in table [l] were used for simulat- 
ing 10 MHz sine wave, 10 MHz square wave and 50 MHz 
sine wave current pumping demonstrating the robustness 
of the simulation. 

From the simulations, we extract an effective channel 
area on the order of 30 nm^ for the pumping (the effective 
channel area is difficult to determine precisely so the values 
are rounded to the nearest lO's of nm^). This was consis- 
tent for both samples SI and S4 shown in figure [l] For 
clarity of presentation, only the data and fits from sample 

51 are shown in figures [2] and |3] Sample S2 and S3 do not 
have a discernible Andreev parameter as it is suppressed 
by their high Ec- This is in general agreement with the 
effective area of 30 nm^ found in earlier work 20 . These 



values are approximately one order of magnitude larger 
than the theoretical channel area of approximately 2 nm^, 
indicating that roughly only one tenth of the effective area 
of the junction is active in agreement with previous results 



31 32 



We have shown that the enhanced current is dependent 
on Ec, pumping amplitude and pumping waveform shape 
but independent of pumping frequency and bias voltage. 
These are all characteristic signatures of Andreev tun- 
neling that can be well accounted for by our theoretical 
model. 

Conclusion. — We have observed the Andreev tun- 
neling process in current pumping with a single electron 
turnstile. This error process can be effectively suppressed 
with high charging energy leading us one step closer to a 
quantum current standard and completing the quantum 
metrological triangle. Andreev reflection as an error pro- 
cess can be fully eliminated with proper choice of Ec as 
demonstrated by these experiments. 
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